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sided estimates for the distribution function F njP (x) of the binomial law 
from the paper published by D. Alters and H. Dinges in 1984. These in- 
equalities are universal (valid for all binomial distribution and all values 
of argument) and exact (namely, the upper bound for F njP (k) is the lower 
bound for F n ^ p (k + 1)). By means of such estimates it is possible to bound 
<^ ! any quantile of the binomial law by 2 subsequent integers. 
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Oh; 

Let X n p be a random variable having the binomial distribution with parameters 

P{X n , p ^k} = J2 c™ P m (i- P ) n - m . 

qq The computation of binomial sums for large n being very tedious, the values of bino- 
£2 mial distribution function are approximated usually by means of the Moivre-Laplace 
^ ! theorem: 
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For example, S.N.Bernstein [2] proved that 




where mo, mi are integers, Zk is the root of equation Zk\/2np(l — p) + ^—^-z\ = 
rrik—np+ak for some ctk G (— §, |) and np(l— p) ^ 62, 5, ^ zq < z\ ^ \j2np{\ — p); 
this formula is valid if max{|mo — np\, \m\ — np\} = O (yjnp{l — p)j . W. Feller [3] 
find a modification of this formula by means of another choice of zq, z\. There are a 
lot of other results on the binomial law, see, e. g. [I]. Relative errors of approximations 
for the tails of binomial distribution function are large due to their superexponential 
decreasing. 

Here we present a new form and give a short full proof of some results due to 
D. Alters and H. Dinges pQ; these authors have used some hints from [5], |EJ. The 
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character of these results is analogous to the Bernstein and Feller theorems, but the 
formulas are explicit and (from the practical viewpoint) constitute almost final solution 
of the large deviation problem for the binomial law. An article pQ remains almost 
unknown (maybe because its presentation is hard to read and the proofs are too long 
and contain nontrivial gaps). Our proof is based on the ideas from pQ. 

Theorem. Let H(x,p) = xln| + (1 — x)ln^| ; sgn(x) = A for x ^ and 
sgn(O) = ; let increasing sequences {C n ^(k)}^ =Q are defined as follows C niP (0) = 
(l-p) n , C niP (n) = l-p n , 

C n , P (k) = $ (sgn (* -p) ^2nH(^p)^j , 1 < k < n. 

Then for every k = 0, 1, . . . , n — 1 and for every p G (0, 1) 

C n , p (k) ^ P{X n , p OX C n , p {k + 1) 

and equalities may happen for k = or k = n — 1 only. 

To demonstrate the accuracy of inequalities (00) 
C n ^ p (k) + C n ^i- P (n — k) = 1; then from 

C nJ) (k) < P{X UjP ^ k}, C n>1 - P (n - k) < P{X n ^ p ^ n 

it follows that 

1 = C n>p (k) + C n ^ p {n - k) < P{X n , p ^k} + P{X n , p ^ k} = 1 + P{X n , p = k). 

In the last inequality the difference between right and left sides is equal to the local 
probability of the binomial law. So, P{X UiP ^ k] — C n:P (k) < P{X UiP = k}. 

The ratio of upper and lower bounds for P{X n ^ p ^ k} in (QQ) may be large if k 
is significantly less than np, but for such k the ratios P{X UiP ^ k + 1}/P{X„^ ^ k] 
are large also. 

In somewhat another form we use the results of [1] in [7] to estimate the partial 
sums of binomial coefficients. 

Proof. Lower bound in (DD) for k = and upper bound for k = n — 1 are exact 
equalities; so these cases will not be considered further. 

We have P{X ni i ^ k] = 0, k < n, and for every integer k G {0, 1, . . . , n — 1} 

k /* 1 J 

P{X»* < £;} = £ C> m (l - ^ m = - / - £ Cn^i 1 ~ z T~ m dz = 

m=0 P m=0 

pi k 

= -n ^2{C^riZ m -\l - zf- m - C^ lZ m {l - z) 71 -™- 1 ) dz = 

JP m=Q 

pi pi 7 k+l(-i _ \n-k-l 

= nC k n _ x \ z k (l- z) n - k - l dz = {k + l)C k+l / { - — -1 dz. 

Jp Jp z 



(i) 



we may note that 



k} = P{X n , p ^ k} 
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Applying the Stirling formula n\ = \/2wn (-J e Sn we find that 
(k + l)C k n +1 = ex P {S k n +1 }> 



(k + l)n n 



2tt(ti - k - 1) (k + - - ' 

where S k+1 = S n - S k+1 - S n - k -i. 

The main step is the proof of the upper bound in (CD) for ^ k ^ n — 2. Denoting 
a = e (0, 1) we have 

n z k+1 (l- z™-*" 1 - ' UZ » f n ( 1_2; M 



= (£) (l^) = exp { ~ n ( a ln 7 + (1 " a) ln fcl) } • 

So, if 5(» = a ln f + (1 - a) ln ±Ef then 

For each a G (0, 1) the function B(z) decreases monotonically from +oo to on (0, a] 
and increases monotonically from to +oo on [a, 1), indeed: 

B(a) = and *(*) = -_ + _ = __. (2) 

It follows that the equation B(z) = ^a 2 (z) has solution a(z) = (a — z)y -^z^z such 
that a(z) decreases monotonically on (0, 1) from +00 to —00. So, 

P{X n , p ^k} = 

- ^ .[ ^ ¥ - ^ ^J? ^ d i • (3) 

Further, let us consider the difference 

5(p) = P{X n , p <: k} - <D (a{p)y/n) , V € [0, 1], 

where $(•) is the standard normal distribution function. IfO^/c^n — 2 then 
5(0) = 5(1) = 0. Now to prove that 5(p) < for all p G (0, 1) it is sufficient to 
show that: a) the function 5(p) is differentiate with respect to p, b) the equation 
5'(p) = has unique root po on (0, 1) and c) 5'(p) < for p G (0, po) and > 
for p G (p , 1). 
In view of 
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Further, it follows from B(p) = \a 2 (p) and (j2j) that B'(p) = = a(p)a'(p), i.e. 

|- $ (o(p) Vn) = <p (o(p)>/n) o'(p) Vn = y> (o(p) Vn) ^fZTp^) = 

We have 



= > ^ r 1 - V^r^fcrJ • (4) 

The multiplier before parenthesis is negative, the first term in parenthesis doesn't 
depend on p. From the formula B(p) = a In + (1 — a) In (^jE^j it follows that 



(p — Cm) 2 / (p — Of) 2 

lim \ I -7—— — - = 0, lim W - NO ^- : =00. 

p|o V (1 - p) 2 #(j9) Pti V C 1 - P? b (p) 

Let us show that the function 2(i-^>fB(p) * s mon °tonicaHy increasing. Its derivative 
equals to 

d (p — a) 2 1 d (p — a) 2 (p — a) 2 d 1 



dp(l—p) 2 B(p) B(p) dp (1 — p) 2 (l—p) 2 dpB(p) 
(p — ck)(1 — a) (p — a) 2 p — a 



1 - pfB(p) (1 - p) 2 p(l - p)B 2 {p) 
2(1 -a)(p- a) ( (p-a) 2 \ 



= 2- 



1 -pfB 2 (p) V 2p(l-a 

The first multiplier in the right hand side changes its sign from — to + at p = a 7 the 
difference in the parenthesis equals to for p = a : and 

— ( B( )- ^ Z a ^ ^ = p ~ a _ p2 ~ a<1 = (P + 1 )(P- a ) 2 > 
5p V 2p(l-a)J Ki-rf 2p 2 (l-a) 2p 2 (l - p)(l - a) > 



for p £ (0, l)\{a}, i. e. the difference — 2 p(i Q ) changes its sign from — to + at 
p = a also. It means that J~ 7iz^g7^ > for all p G (0, l)\{a}. 

Consequently, the difference eS " l -\/j~~ \j 2{i^yW ~p) on IP' 1] decreases monoton- 
ically from e Sn to —00 and is equal to for some p = pq. We see that 5(p) < 
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for all p G (0, 1), i. e. for all n ^ 1, ^ k ^ n - 2, p G (0, 1) 

P{X n , p ^ £;} < $ (o(p)Vn) = 

= $ (sgn (*±1 - p) ^2ntf = C n , p (fc + 1). 

Now we prove the lower bound in (CD) for k G {1, . . . , n — 1}. As 

P{X n?p ^ £;} = P{X M _ P ^n-k} = 1- P{X M _ P ^ n - k - 1} 
and n-fc-1 E {0, . . . , n — 2} we have (using upper bound just proved) 

P{X n , p ^ k} > 1 - <D (sgn (S=* - (1 -p)) ^2nH , 1 - ^ = 
= $ (sgn (J - p) y/2nH(^ P )^ = C n , p (k). 

The theorem is proved. 

Remark. Inequalities (OQ) may be sharpened by means of nontrivial upper bounds 
of the function 
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For concrete values of parameters this integral may be estimated numerically. 

To obtain the analytic estimates we may note that the integrand is the prod- 
uct of continuous positive function f(z) = - (p (a(z) y/n) \Jn and monotonically de- 
creasing function g(z) = e S " i+1 — \J 2{i^zfB {z)^ an< ^ that Jo f( z )9( z )^ z = 0- So, 
if po G (0, 1) is such that g(po) = 0, then the value of the integral may be bounded 
from below: for < p ^ po we have 



rp pit 

/ f( z )g( z )dz ^ sup g(u) / f(z)dz, 

'0 0<u^p JO 



and for po $C p < 1 
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f(z)g(z)dz = - / f(z)g(z)dz^ sup / /(z)^. 



These estimates along with (j3j) and (OQ) give for p ^ po 



S(p) 



#(y) dy ^ 



= — sup 

*G(p,l) 



fc+l / 



sup 

zG(p,l) 



.2 - a 



1-a V 2 ( X - z) 2 B{z) 



f 1 

/ -V (a(y)Vn) Vndy = 

z y 



2a(l - £) 2 £(*) 



^ — sup 

«€(p,l) 



1 - e 



i j {l-a){z-a) 2 
2a{l - z) 2 B{z) 



Cn 7 z(k). 
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Analogously, for p < po 



Jo 



= — sup 

ze(o, P ) 



jfc+i / ex 



z — a) 



1-a \ 2(1 -z) 2 B( 



= — sup 

2G(0jj) 



^ — sup 

zG(0,p) 



1 - e 



-s 



fc+i (1 — a)(z — a) 2 



2a{l - z) 2 B{z) 



f z 1 

/ -V (a{y)Vn) Vndy = 

'o y 



P{X n , z > fc + 1} < 



1 - e 



-S* 



l / (1 — a)(z — a) 2 



2a{l-z) 2 B{z) 



'l-C^k + l)). 



(6) 
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Or 

Let X n p be a random variable having the binomial distribution with parameters 
g; (n,p): 

&. P{X n , p ^k} = C™p m (l-p) n - m . 

qq Computation of binomial sums for large n being very tedious, the values of binomial 
£2 ; distribution function are approximated usually by means of the Moivre-Laplace theo- 
^ ! rem: 

lim P \X np ^ np + x\Jnp(l — p) \ = &(x) = / ip(u) du, tp(u) = e~ u ^ 2 . 

n-Hx> I J J ' V27T 
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For example, S.N.Bernstein [2] proved that 



E 



mi— 1 , , 

c k nP k (i- P y 

=m 




where mo, mi are integers, Zk is the root of equation Zk\/2np(l — p) + 



rrik—np+ak for some ctk G (— §, |) and np(l— p) ^ 62, 5, ^ zq < z\ ^ \j2np{\ — p)\ 
this formula is valid if max{|mo — np\, \m\ — np\} = O (yjnp{l — p)j . W. Feller [3] 
modified this formula by a different choice of zq, Z\. There are many other results 
on the binomial laws, see, e. g. [I]. Relative errors of approximations for the tails of 
binomial distribution function are large due to their superexponential decreasing. 

Here we present a new form and give a short full proof of some results due to 
D. Alters and H. Dinges [1]; these authors have used some hints from [5], |5]. The 
character of these results is analogous to the Bernstein and Feller theorems, but the 
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formulas are explicit and (from the practical viewpoint) constitute almost final solu- 
tion of the large deviation problem for the binomial law. Results [T] remains almost 
unknown (maybe because this article is difficult for reading, the proofs are too long 
and contain nontrivial gaps). Our proof is based on the ideas borrowed from [T]. 

Theorem. Let H(x,p) = xln| + (1 — x)hij^ } sgn(x) = r| for i^O and 
sgn(O) = ; let {Cn 7 p(k)}^ be increasing sequences defined as follows: C njP (0) = 
(l-p) B , C n , p (n) = l-p n 7 

C n , P (k) = $ (sgn (* -p) ftnHfop)^ , 1 < k < n. 

Then for every k = 0, 1, . . . , n — 1 and for every p G (0, 1) 

C n , p (k) ^ P{X n , p OK C n , p (k + 1), (1) 

and equalities may happen for k = or k = n — 1 only. 

To demonstrate the accuracy of inequalities (0Q) we observe that 
,i- p {n - k) = 1. Then it follows from 

C nj> (k) < P{X UjP ^ k}, C n ,i_p(n - k) < P{X nA ^ p <^n-k} = P{X UjP ^ k} 

that 

1 = C njP {k) + C n ^ p {n -k)< P{X n , p ^k} + P{X n , p ^k} = l + P{X n , p = k}. 

In the last inequality the difference between the right and left sides is equal to the 
local probability of the binomial law. So, P{X UjP ^ k} — C n ^ p {k) < P{X UjP = k}. 

The ratio of the upper and lower bounds for P{X n p ^ k] in (0Q) may be large if k 
is significantly less than up, but for such k the ratios P{X UiP ^ k + l}/P{X n ^ p ^ k} 
are large also. 

In a different form we use the results of [I] in [7] to estimate the partial sums of 
binomial coefficients. 

Proof. The lower bound in (DQ) for k = and the upper bound for k = n — 1 are 
exact equalities; so these cases will not be considered further. 

We have P{X n i ^ k] = 0, k < n, and for every integer k G {0, 1, . . . , n — 1} 

k /* 1 j 

P{X n , p O} = £ C™ P m (l - P ) n - m = - / ^ E Cn^i 1 ~ dz = 

m=0 P m=0 

„1 k 

= -n ^2{C^riZ m -\l - zf- m - C^ lZ m {l - z) 71 -™- 1 ) dz = 

JP m=Q 

= nC k n _ x \ z k (l- z) n - k - l dz = {k + l)C k+l / { - — -1 dz. 

Jp Jp z 
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Applying the Stirling formula n\ = \/2wn (-J e Sn we find that 
(k + l)C k n +1 = exp{S* +1 }> 



(k + l)n n 



2tt(ti - k - 1) (k + l) fc+1 (n - - 1) 



n— fc— 1 



where = S n - S k +i - S n -k-i- 

The main step is the proof of the upper bound in (CD) for ^ k ^ n — 2. Denoting 
a = *±! € (0, 1) W e have 

T7 / \ / /-I \ \ 71— fc— 1 



(A; + l) k+1 (n - k - v 7 \A: + iy V^-^" 1 

= (aJ — =eX p -n aln- + (l-a)ln — 



So, if = a In f + (1 - a) In then 

y 27r(l - a) J p z 

For each a G (0, 1) the function B(z) monotonically decreases from +oo to on (0, a] 
and monotonically increases from to +oo on [a, 1). Indeed, 

B(a) = and B'(z) = + = —. r . 2 

z 1 — z z(l — Z) 

Hence it follows that the function a(z) = (a — z)J solves the equation 
B(z) = \a 2 (z) and decreases monotonically on (0, 1) from +oo to — oo. So, 

P{X n , p ^k} = 

Further, let us consider the difference 

5(p) = V{X niP <: k] - <D (a(p)v^) , P € [0, 1], 

where $(•) is the standard normal distribution function. IfO^/c^n — 2 then 
5(0) = 5(1) = 0. Now to prove that 5(p) < for all p G (0, 1) it is sufficient to 
show that: A) the function 5(p) is differentiate with respect to p, B) the equation 
5'(p) = has a unique root po on (0, 1) and c) 5'(p) < for p G (0, po) and 5'(p) > 
for p G (po, !)• 
In view of (ED 



d ^ r ^ r , , cfe+i / 1 / , . . — x 
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Further, it follows from B(p) = \a 2 (p) and (j2j) that B'(p) = = a(p)a'(p), i.e. 

a'(p) = pi i P S p )a( P ) > and therefore 

|- $ (o(p) Vn) = <p (o(p)>/n) o'(p) Vn = (p {a{p)^) ^^ p ^ = 

, r , yjn / (p - a) 2 

We have 



The multiplier before the parenthesis is negative, the first term in the parenthesis 
doesn't depend on p. From the equality B(p) = a In + (1 — a) In f^E^l ^ follows 
that 

(p— a) 2 (p ~ a ) 2 
hm ^ / — - . . = 0, lim h / NO - . - = oo. 

pid V (i - p) 2 b(p) P n y (l - p) 2 £(p) 

Let us show that the function 2 (i P p )^B(p) m onotonically increases with p. Its derivative 
equals to 

d (p — a) 2 1 d (p — a) 2 (p — a) 2 d 1 



dp(l—p) 2 B(p) B(p)dp(l—p) 2 (l—p) 2 dpB(p) 
(p — a)(l — a) (p — a) 2 p — a 



= 2- 



(1 - pfB{p) (1 - p) 2 p(l - p)B\p) 
2(1 - a) (p-ot) ^ (p-a) 2 \ 



[\-pfB\p) V 2^(1- a) / ' 

The first multiplier in the right hand side changes its sign from — to + at p = a, the 
difference in the parenthesis equals to for p = a, and 

— ( b( ) - ^ - a ^ 2 ^ = ^~ Q _ ^ 2 ~ ^ 2 = b + «) 2 

dp V 2p(l - a)/ p(l - p) 2p 2 {l - a) 2p 2 {l - p){l - a) > 



for p G (0, l)\{a}, i. e. the difference B(p) — 2 p(i a ) changes its sign from — to + at 

d (p-a) 2 
' dp(l~p) 2 B(p) 



p = a also. Thus, Ir yAAy > for all p G (0, l)\{a}. 



Consequently, the difference e 5 ™ y/j^ ~ y 2(1 on t^' ^ * s monotonically 

decreasing from e^«y / ^^ to —00 and is equal to for some p = po- We see that 
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6{p) < for all p 6 (0, 1), i. e. for all n ^ 1, ^ k ^ n - 2, p G (0, 1) 

P{X n , p ^ £;} < $ (o(p)Vn) = 

= $ (sgn (*±1 - p) ^2ntf (^.p)) = C n , p (fc + 1). 

Now we prove the lower bound in (CD) for k G {1, . . . , n — 1}. As 

P{X n?p ^ £;} = P{X n ,i_ p > n — fc} = 1 — P{X n ,i_ p ^ n - A; - 1} 
and n — k — 1 G {0, . . . , n — 2}, we have (using the upper bound just proved) 

P{X n , p ^ k} > 1 - <D (sgn (S=* - (1 -p)) (i^*,l-p)^ = 

= <D (sgn (J - p) y/2nH(^ P )^ = C n , p (k). 

The theorem is proved. 

Remark. Inequalities (CD) may be sharpened if nontrivial upper bounds of the 
function 



'0 



' J „ (a(z)VS) K J+ '\/rS " J 2(1 ( : K P e (0. !)■ 



is known. For concrete values of parameters this integral may be estimated numerically. 

In the general case we observe that the integrand is the product of a continu- 
ous positive function f(z) = - ip (a(z)y/n) y/n and a monotonically decreasing func- 
tion g(z) = e s " +1 - \l 2{ i_~yb {z ) , and that J* f(z)g(z)dz = 0. So, if p G (0, 1) 
is such that g(po) = then the value of the integral may be bounded from below: for 
< p ^ po we have 



'o 



f(z)g(z)dz ^ sup g(u) / f(z)dz, 

0<u^p JQ 



and for po ^ p < 1 
rP 



f(z)g(z)dz = - / f(z)g(z)dz^ sup / f(z)dz. 

'0 ./jD P<M<1 Ju 



These estimates along with (j3j) and (OQ) give for p ^ po 



S(p) 



#(y) dy ^ 



= — sup 

*G(p,l) 



fc+l / 



sup 

zG(p,l) 



.2 - a 



1-a V 2 ( X - z) 2 B{z) 



f 1 

/ -V (a(y)Vn) Vndy = 

z y 



2a(l - £) 2 £(*) 



^ — sup 

«€(p,l) 



1 - e 



i j {l-a){z-a) 2 
2a{l - z) 2 B{z) 



Cn 7 z(k). 



(5) 



Analogously, for p < po 



Jo 



= — sup 

ze(o, P ) 



jfc+i / ex 



z — a) 



1-a \ 2(1 -z) 2 B( 



= — sup 

2G(0jj) 



^ — sup 

zG(0,p) 



1 - e 



-s 



fc+i (1 — a)(z — a) 2 



2a{l - z) 2 B{z) 



f z 1 

/ -V (a{y)Vn) Vndy = 

'o y 



P{X n , z > fc + 1} < 



1 - e 



-S* 



l / (1 — a)(z — a) 2 



2a{l-z) 2 B{z) 



'l-C^k + l)). 



(6) 
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